Water Transport through Carbon Nanotubes with the Radial Breathing Mode 
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Molecular dynamics simulations are performed to investigate the water permeation across the 
single-walled carbon nanotube (SWCNT) with the radial breathing mode (RBM) vibration. It is 
found that the RBM can play a significant role in breaking hydrogen bonds of the water chain, 
accordingly increasing the net flux dramatically, and reducing drastically the average number of 
water molecules inside the tube with the frequency ranging from 5000 to IIOOO GHz, while far away 
from this frequency region the transport properties of water molecules are almost unaffected by the 
RBM. This phenomenon can be understood as the resonant response of the water molecule chain 
to the RBM. Our findings are expected to be helpful for the design of high-flux nanochannels and 
the understanding of biological activities especially the water channelling. 

PACS numbers: 31. 15. at, 31.15.xv, 61.20.Ja 



Introduction — The transport properties of water 
molecules across nanochannels have been extensively in- 
vestigated in the past [l|-l7[. In recent years, various 
designs of nanochannels concerning temperature gradi- 
ents [8[, charge modification [9|, coulomb dragging [10[, 
and static electric fields |ll[ have been proposed to study 
abundant novel properties of water permeation. These 
progresses may have significant technical implications for 
the design of novel nanofluidic devices or machines, such 
as desalination of seawater 12 1, molecular sieves [13[, 
molecular water pump |14l | , and so on. Furthermore, one 
important reason why research on nanochannel transport 
properties has continued to be of interest is that some pri- 
mary characters of water confined in simple nanochannels 
are similar to that of complex biological channels. Typi- 
cally, water across both the biological membrane proteins 
(for instance, aquaporin-1 (AQPl) and Glpf) [Ig] and 
the nanochannels [l6[ can form single-file water chains of 
stable dipolc orientation. The vapor-liquid transition has 
also been observed in mechanoscnsitive membrane chan- 
nels 17 1 and other hydrophobic nanopores [l8[ . Another 



resemblance existing in the nanochannels and biological 
channels is that they can share both wavelike density 
distribution and wet-dry transition resulting from con- 
finement [li,!!^. 

It has been long well-known that the wave mechanics 
well interprets the truth of nature. Microscopic objects 
such as biological microtubules are actually exposed in 
the environmental vibration modes of all possible micro- 
oscillators entangled with intrinsic modes. It is thus of 
broad interest to study the molecular permeations with 
the vibrational coherence and/or interference. Indeed, 
theoretical and experimental studies indicate that the 
single-walled carbon nanotubes (SWCNTs) have various 
collective vibrational modes |2lM2J|. Collective motion, 



different from the internal fast processes, is usually re- 
garded as the fifth dimensional motion, whereas it re- 
mains almost untouched in the study of water perme- 
ation through SWCNTs. The radial breathing mode 
(RBM) [25|, being the simplest collective vibration, is 
one of the most important modes of SWCNTs, and it can 
be excited easily by radio wave techniques or constantly 
by environmental interferences. In this study, wc thus 
restrain ourselves to the RBM. It is known that molecu- 
lar dynamics (MD) simulations have been widely used to 
investigate the transport pro perties of water molecules 
across SWCNTs |3, S H lli [1^ . In this work, we use 
MD simulations to explore the dynamic behavior of water 
permeation through SWCNTs with the RBM. We have 
observed that the RBM of the SWCNT can greatly af- 
fect the water permeation properties through a resonant 



Computational Methods — The simulation framework 
is illustrated in Fig. [TJ An uncapped (6,6) armchair 
SWCNT with a length of 1.34 nm and a diameter of 
0.81 nm is embedded in two graphite sheets along the z 
direction. The distance between the bottom end of the 

o 

SWCNT and the graphite sheet is 2 A. Periodic bound- 
ary conditions are applied in all directions. The elec- 
trostatic calculation, carried out at a constant volume 
with box size dimensions of L^; = 3.5 nm, Ly = 3.5 nm, 
Lz = 6.3 nm and constant temperature (300 K) achieved 
by Langevin dynamics, is done using the PME 127| sum- 
mation method with a cutoff for real space of 1.2 nm. The 
cutoff for the van der Waals interaction is also 1.2 nm. 
With the MD program NAMD2 ^, all simulations are 
performed using the CHARMM27 force field [Uls^ and 
the TIP3P water model [3l[ . The time step of Ifs is used, 
and data arc collected every 0.5 ps. To obtain a directed 
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flow, an additional acceleration of 0.01 nmps~ , equiva- 
lent to a pressure difference about 20 Mpa between two 
ends of the SWCN'R is applied to each water molecule 
along -l-z direction [3|, [SOJ . 




FIG. 1: (Color online) Snapshot of the simulation system. 
An uncapped (6,6) armchair SWCNT combined with two 
graphite sheets is solvated in a water box of 3.5x3.5x6.3nm^ 
with 1852 molecules. There is a vacuum between two sheets. 
The axis of SWCNT is parallel to the z-axis. 



The RBM is simply simulated with a coherent radial 
motion obeying the following relation 

r{x, y, t) = ro(x-o, yo, 0) + Acos{ujt + ip), (1) 

where rg is the equilibrium position of each carbon atom, 
and A is the amplitude. In the calculation, we take 
ip = 7r/2 and A to be 5 percent of the intersection ra- 
dius. This percentage is significantly larger than the one 
in the literature J32M35| . and we choose it for the fol- 
lowing considerations. Though this work focuses on the 
RBM effect, one may be more interested in the energy 
transfer or molecule transport process of the microsys- 
tems induced by the more general and complex collective 
vibrations that involve the RBM and usually have much 
larger amplitudes [2l|, 



36, 371. On the other hand, the 



vibration amplitude of a biological microtubule can also 
be much larger 40|, |41|. Specifically, we have tested that 
the (6,6) SWCNT is rather fiexible and it is not very dif- 
ficult to reach the vibration amplitude as much as 5 % by 
appropriate external induction 38|, |39| . Without losing 
the generality, we thus choose a moderate amplitude 5% 
to observe clear effects of collective vibrations and pro- 
vide useful insights for the simulations of the biological 
microsystems. Moreover, because the filling of the water 
just affects the RBM frequency slightly [421, we neglect 
it in our estimation. 

Results and Discussion — The simulation is performed 
for the SWCNT system with different RBM frequencies. 
For each frequency, the simulation time is 110 ns. Un- 
less otherwise indicated, the data of the last 100 ns is 
used for analysis. For clarity, we define the net flux 
to be the difference between the water molecular num- 
ber per nanosecond across the SWCNT from one end to 
the other [3|, |l9[ . The average number [3|, ll6| of water 
molecules inside SWCNT is denoted with the symbol N. 

Fig. [5] shows the net flux and N of water molecules 
across the SWCNT as a function of RBM frequency. 
It is observed that the net flux and N do not change 
much at most frequencies and are about 20 ns^^ and 5, 



respectively. However, as shown in Fig. [2] (inset), the 
net fluxes undergo a clear enhancement in the frequency 
range from about 2000 GHz to 14000 GHz, and two peak 
values around 4000 GHz and 12000 GHz are 63 ns'^ and 
55 ns~^, respectively, being approximately three times 
those beyond this frequency range. As shown in Fig. [2j 
we observe that the N is around 1 in the frequency range 
of 5000-10000 GHz, and its minimum is about 0.84 at 
6000 GHz. Water molecules inside the SWCNT are ac- 
tually absent within 5000-10000 GHz in a majority of 
time. As an example, we have counted up the temporal 
distribution of the average number of water molecules 
inside the SWCNT in the 100 ns at 10^ GHz. We find 
that the empty state, namely, no water molecule inside 
the tube, accounts for 46.1% of total time, and the oc- 
cupancy with N = 1,2,3,4, and 5 accounts for 22.5%, 
14.5%, 9.9%, 5.5%, and 1.5%, respectively. Meanwhile, 
we have counted the hydrogen bond number of water 
molecules inside the SWCNT in this frequency range, ac- 
cording to the criterion of the oxygen distance less than 

3.5 A and hydrogen-bond angle < 30° [19[, and consis- 
tently found that it is almost always zero. Clearly, with 
hydrogen bonds of water molecules being almost com- 
pletely broken, the water molecules move individually in 
the SWCNT with the RBM in this frequency range, in 
contrast to the existence of the single-file water chain 
at other frequencies. Consequently, the sharp enhance- 
ment of the flux follows from the full fracture of hydrogen 
bonds and fast shuttling of water molecules from one end 
to the other. 
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FIG. 2: (Color online) The net flux and average number A'' 
of water molecules inside the tube as a function of RBM fre- 
quency. Given in the inset is an image amplification within 
the frequency interval 0-16000 GHz. 

A coincident breaking of the hydrogen bonds can be 
realized by a perturbation of the form e*"* in suffi- 
ciently long time span, while the understanding of a thor- 
ough fracture can resort to a classical harmonic oscillator 



model. Supposing the intrinsic frequency of the single- 
file water chain is determined by the binding energy of 
the hydrogen bond, Eh, we can estimate the classical 



resonant frequency using the relation Eh 



2^2/2 



where m is the mass of water molecule, and A is the 
amplitude of allowed radial motion. Using the value of 
Eh ~ 16 kcal/mol [16[ and the amplitude about 0.3-0.9 

o 

A which are half-peak width values of the Gaussian dis- 
tribution observed in the simulation, the classical res- 
onant frequency of water chain is estimated to lie be- 
tween 4800 GHz and 14500 GHz which is well situated 
in the frequency range for large fluxes in our simulations. 
It is noteworthy that the theoretical value of inherent 
frequency of the short (6,6) SWCNT is about 10* GHz 
[nisi, Sill. With the REM frequency away from the 
resonance, segments of the water chain instead of indi- 
vidual molecules are produced and hence the occupancy 
of water molecules increases. A double-peak structure of 
the water flux thus appears due to the fast shuttling of 
the broken segments. With the REM frequency far away 
from the resonance, the flux is almost unaffected by the 
REM because there is just coincident breaking of hydro- 
gen bonds. Experimentally, this drastic resonant effects 
on the water permeation can be verified with radio wave 
techniques. 

To better understand the process of the hydrogen bond 
fracture, a dynamic evolution of the average kinetic en- 
ergy of water molecules held in the tube is shown in 
Fig. |3l It is interesting to see that the significant fluc- 
tuation at / = 10* GHz starts to form actually at a 
very early time and then the shift to a strong oscillation 
occurs after a relatively long latent period. While this 
phenomenon does not take place at other frequencies as 
shown in Fig. [3j it can be understood as a resonant re- 
sponse to the REM of the SWGNT. We note that this 
resonance accompanied with the H-bond rupture needs 
an energy accumulation process with the time span de- 
pending on the REM amplitude. 

Shown in Fig. 2] is the distance of the oxygen atom 
of the water molecule inside the SWCNT from the tube 
axis evolving with the time. Similarly, we see again that 
the oscillation is forced to resonate with the REM of the 
SWCNT at 10* GHz, while the coherent oscillations do 
not occur at other frequencies. Posterior to the fracture 
of hydrogen bonds, the radial oscillation is dictated by 
the van der Waals interaction between the carbon atoms 
of the SWCNT and water molecules. It is known that 
the motion of water molecules is driven by the potential 
field. We thus calculate the average electrostatic, van 
der Waals and total interaction potentials: Pe, Py and 
Pt ^ Pe + Py per water molecule at different position 
along z axis inside the tube where the water molecule 
interacts with the SWCNT, two graphite sheets, and all 
ambient water molecules. As shown in Fig. [51 the av- 
erage interaction potentials are almost the same at all 
frequencies but / = 10* GHz. This explains why the net 
flux and N do not have significant differences at most 
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FIG. 3: (Color online) The kinetic energy for each water 
molecule inside the SWCNT with respect to time. Results 
in an earlier period are presented in the lower panel where 
the colored curves are displaced downwards for clarity. 
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FIG. 4: (Color online) The distance of the oxygen atom of the 
water molecule from the z-axis evolving with time for various 
RBM frequencies. 



frequencies. In fact, we can continuously observe sta- 
ble one^imensionally ordered water chain like previous 
16l Il9j using molecular visualization program 
at these frequencies. However, at 10* GHz, 
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Fig. [5] (a) shows that the Pe is close to kcal/mol, which 
means that hydrogen bonds of water molecules inside 
the SWCNT are almost completely deprived. Fig. [5] (b) 
shows that the Py is repulsive for single water molecules, 
revealing the fact that the single water molecules are 
forced to move closely to the tube wall. The Pj, as shown 
in Fig. \5lc), demonstrates a similar character featured 



by Pf, and P^. Due to the repulsion provided by the 
RBM resonance, it is difficult for water molecules to en- 
ter into and keep in the SWCNT, resulting in a dramatic 
reduction of N , whereas this can indeed provide an ef- 
ficient transport of water molecules accompanied by a 
large net flux. We mention that the peak values of the 
net flux do not coincide with the minimum values of N 
as discussed above. On the other hand, wc note that 
the water molecules can be as far as 2 Angstrom from 
the z axis when the resonance happens (see Fig. 0]). In 
such cases, it is necessary to explore the validity of the 
force-field model. With comparison to the total ener gy o f 
the snapshots obtained from the ab-initio approach 



we find that the energy acquirement of water molecules 
is clearly overestimated by the MD simulation. By re- 
ducing the radius parameters in the MD simulation, we 
may roughly reproduce the ab-initio results. In doing so, 
the phenomena of resonance remain unchanged, while the 
flux is just moderately enhanced. Wc observe that the H- 
bond fracture due to the resonant mechanism are mainly 
relevant to the relatively long-range van der waals inter- 
action in the vicinity of the equilibrium position, while 
the short-range repulsive interaction plays a more effi- 
cient role posterior to the H-bond fracture. Therefore, 
the main conclusions that we have drawn are not essen- 
tially affected by the deficiency of the MD method. 
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FIG. 5: (Color online) The average electrostatic (a), van der 
Waals (b), and total (c) potentials per water molecule that 
interacts with the SWCNT, two graphite sheets, and all am- 
bient water molecules as a function of z-coordinates. 

It is worth discussing the dependence of the flux on the 
RBM amplitude. The RBM effect on the flux decreases 
with the decrease of the RBM amplitude. For instance, 
for an amplitude of 3% of the intersection radius, the 
flux at / = 10"' GHz reduces to be about 40 ns''^ which 
is just twice the value at frequencies without the RBM 



influence. The flux enhancement can further be reduced 
till to disappearance at a less RBM amplitude. With 
the present setup, the smallest RBM amplitude is about 
2.7%, while it can reduce moderately with decreasing the 
imposed pressure gradient. 

At last, it is significant to explore the RBM effects on 
the water transport in SWCNTs with generalized setups 
that may bridge the connection with the experiments us- 
ing the larger diameters and longer lengths. We have 
chosen the SWCNT setups with the diameter ranging 
from 1.09 to 1.36 nm and the length from 1.34 to 10 
nm. Likewise, we have observed that the rupture of the 
H-bonds of various hinged water chains can occur in all 
these SWCNTs with the RBM at resonant frequencies, 
resulting in a noticeable difference in water fluxes. More- 
over, we have studied the RBM effects on water trans- 
port properties at much smaller pressure gradients. It 
is interesting to observe that the water chain fracture is 
more remarkable for a small pressure gradient, e.g., 15 
atm, owing to a much longer stay of water molecules in 
the tube. These investigations indicate that the RBM 
resonance is rather universal. 

Conclusions — In summary, we have demonstrated ef- 
fects of the resonant mechanism, responding to the RBM 
of the SWCNT, on the transport properties and dynami- 
cal behaviors of water molecules across the SWCNT. The 
resonance of the single-file water chain with the RBM of 
the SWCNT, characterized by the full fracture of hy- 
drogen bonds, results in a dramatic decrease of water 
molecules inside the tube and a sharp enhancement of the 
net flux. It is also found that away from the resonance on 
the both sides roughly in a frequency range 2000-14000 
GHz the RBM of the SWCNT can create partial breaks 
of hydrogen bonds and accordingly dramatic changes in 
transport properties of water molecules, whereas the lat- 
ter are almost unaffected beyond this range. In this fre- 
quency range, a double-peak structure of the net flux 
forms with the vast enhancement of the net flux, and the 
peak values around 4000 GHz and 12000 GHz are three 
times those unaffected by the RBM. Our findings may 
have some possible implications to biological activities 
and potential technical applications in the development 
of nanotube-based nanofluidic devices. 
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